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Abstract 

Q"^ \ We investigate dynamics of the chiral transition in expanding quark-antiquark 

plasma produced in an ultra-relativistic heavy ion colhsion. The chiral symmetry 
break-down and dynamical generation of the constituent quark mass are studied within 
(^ \ the linear sigma model and Nambu-Jona-Lasinio model. Time dependence of the 

r^ ■ quark and antiquark densities is obtained from the scaling solution of the relativistic 

Vlasov equation. Fast initial growth and strong oscillations of the constituent quark 
mass are found in the linear sigma model as well as in the NJL model, when derivative 
terms are taken into account. 
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Introduction. - It is commonly believed that colour deconfinement and chiral symmetry 
restoration take place at early stages of ultra-relativistic heavy-ion collisions. At interme- 
diate stages of the reaction the quark-gluon plasma may be formed and evolve through the 
states close to thermodynamical equilibrium. However, at later stages of the expansion the 
transition to the hadronic phase with broken chiral symmetry should take place. The break- 
down of chiral symmetry will possibly lead to such interesting phenomena as formation of 
disoriented chiral condensates (DCCs) and classical pion fields as well as clustering of quarks 
and antiquarks. These phenomena were studied recently in many publications [1-14], using 
QCD motivated effective models, such as the linear and non-linear sigma models and the 
Nambu-Jona-Lasinio (NJL) model. Of course, these models have some significant shortcom- 
ings, e.g. they do not possess colour confinement and the NJL-model is non-renormalizable. 
The key point is, however, that these models obey the same chiral symmetry as the QCD 
Lagrangian. 

In most applications of the sigma model the quark degrees of freedom are disregarded (see 
e.g. [2-9,14]). The inclusion of quarks [|10|, ^ makes it possible to study the hadronization 



process, in particular, the dynamical generation of the constituent quark mass. In this letter, 
we consider the late stages of the plasma evolution, such that collisions between quarks and 
antiquarks are not frequent enough to maintain thermodynamical equilibrium. We assume 
that the expansion is spherical at this stage. We use the linear sigma model and the NJL 
model to describe the interaction of quarks with background chiral fields. 

Linear sigma model.- The Lagrangian density of the linear sigma model with quark 
degrees of freedom reads 

C = q[^^^^^, - g{a + i^^f ■ 7f)]g + ^{d.cxd^'cx + d.nd^'ff) - U{a, tt), (1) 

where 

U{a, 7f) = ^{a^ + n^- v^ - Ha (2) 

is the Mexican Hat potential. Here q is the light quark field q = {u,d). The scalar field a 
and the pion field vf = (tti, 7r2, 773) form together a chiral field $ = (cr, tt). This Lagrangian 
is approximately invariant under chiral SUlC^) (8) SUr{2) transformations if the explicit 
symmetry breaking term Ha is small. The parameters of the Lagrangian are chosen so that 
the chiral symmetry is spontaneously broken in the vacuum and the expectation values of 
the meson fields are (a) = Jt^ and (vr) = 0, where /^r = 93 MeV is the pion decay constant. 
The constant H is fixed by the PCAC relation that gives H = firfn"^, where 171^^ = 138 MeV 

2 

is the pion mass. Then one finds v^ = /^ — ^. The A^ is determined by the sigma mass, 
m^'^ = IX'^fl + m^, which we set to 600 MeV yielding A^ ~ 20. The couphng constant g is 
fixed by the requirement that the constituent quark mass in vacuum, rrivac = gf-K, is about 
1/3 of the nucleon mass, that gives g ~ 3.3. With these parameters a chiral phase transition 
is predicted at T, ^ 132 MeV Q. 

Below we adopt the mean field approximation, considering a and vf as classical fields. 
The variation of C with respect to a and vf yields the equations of motion 

d^d^aix) + \^[a\x) + 7f2(a;) - v'']a{x) -H = -gps{x), 

d^d^^{x) + \^[a'^{x) + 7f2(x) - fVW = -gpp{x), (3) 

where ps = (qq) and Pp = iiq'j^rq) are respectively the scalar and pseudoscalar densities 
of valence quarks. They are determined by the interaction of quarks with meson fields. As 
shown in [|10] these densities can be expressed as Ps{x) = ga{x)a{x) and Pp = g7i{x)a{x), 
with scaling function 



a[x] — '' 



Vq f d^p 



(27r) 



p2 _|_ jji2(^x) 



[nq{x,p)+ng{x,p)] (4) 



Here z/g = 2x2x3 = 12is the degeneracy factor of quarks and m{x) is the constituent 
quark mass which is determined self-consistently through the meson fields, 

m\x)=g^[a\x) + n\x)]. (5) 

In eq. (^ n^ and rig are the quark and antiquark occupation numbers which in thermal 
equilibrium are given by the Fermi-Dirac distributions. 

Nam,hu-Jona-Lasinio model.- Now we turn to the Nambu-Jona-Lasinio model which is 
widely used now for describing hadron properties and the chiral transition. The Lagrangian 
as introduced by Nambu and Jona-Lasinio |ff5l in 1961 is given by: 



G 

C = q(^iryt^d^-mo)q+ — {qqf + {qil^rqf , (6) 



where q stands for the quark field and uiq is the small current quark mass. At vanishing mo, 
the NJL Lagrangian is invariant under the chiral SUl{2) ® SUii{2) transformations. The 
coupling constant G has dimension (energy)"^, that is why the theory is non-renormalizable. 
Therefore, a cutoff momentum A is introduced to regularize divergent integrals. It defines 
an upper energy limit for this effective theory. Free parameters of the model are fixed 
to reproduce correctly the vacuum values of the pion decay constant, pion mass and the 
constituent quark mass. When the current quark mass mo is set to zero, the remaining 
parameters are G = 11 GeV~^ and A = 653 MeV [O]. With these parameters the chiral 



transition occurs at critical temperature T^ = 190 MeV \W\ which is significantly higher 
than in the sigma model. 

In the mean-field approximation the Lagrangian (0) is represented in a linearized form 



17| where the constituent quark mass is expressed as m(x) = ms{x) + i'y^r ■ mp{x). By 
definition 

ms{x) = -G{{-q{x)q{x))), mp{x) = -G {{iq{x)'j5Tq{x))) , (7) 

where ((...)) means averaging over the exact many-body state including the Dirac sea. The 
right hand sides of these expressions can be expressed through the quark Green's function 
S{x,y) which is determined self-consistently with the exact mass rh{x). In the case of slow 
varying mass the scalar and pseudoscalar densities can be parametrized in terms of quark 
and antiquark occupation numbers as Psix) = ms{x)a{x) and p^(x) = rhp{x)a{x) where 

Q(^) = 7T^/ /, ^ , =K(a^,p) + r^g(x,p)-l]. (8) 

{27iy J Jp-i + m^ix) 



Here the term with —1 corresponds to the contribution from the Dirac sea and the rest 
comes from valence quarks and antiquarks (compare with eq. (^)) 

In the case of rapidly changing m{x) one should solve the exact consistency relations in 
terms of S{x, y). For this purpose the derivative expansion method was proposed by Eguchi 
and Sugawara [|1^. It consists in expanding S{x, y) around the Green's function So{x, y) for 
a homogeneous state and retaining only divergent terms in the momentum-space integrals. 



In contrast to |T^, where only the vacuum contribution was considered, here we expand 
around the state characterized by the scalar and pseudoscalar densities Ps{x) and Pp{x). 
Then, with minor modifications, one can repeat all steps of the derivation and arrive at the 
modified gap-equation: 

d^d^rh{x) - 2G^[pI{x) + fp{x)]rh{x) + 2\m{x)\^m{x) = 0. (9) 

Here |m(a;)p = m?{x) = rriKx) + mp(x) and p1{x) + p Jx) = m^(x)a^(x), where a^(x) is 
given by eq. (|). In the limit of a homogeneous system {d^ = 0) this equation is reduced to 
G^a^ = 1 which is equivalent to the famous gap equation used in most applications of the 
NJL model. Eq. (^ can be easily generalized to include the explicit symmetry breaking 
terms. By expressing mo in terms of m^ they can be represented in the form analogous to 
eq. ®. 

One can easily see that in vacuum eq. (^) is equivalent to the linear sigma model 
discussed above. Indeed, in this case Pp = and the gap equation has the form rriyac = —Gps, 
where rriyac is the constituent quark mass in vacuum, which is written above as rriyac = gfn- 
Then eq. (H) coincides with eqs. (H) if one chooses A^ = 2g'^. (In fact this is very close to our 
choice of parameters) . With this choice m„ = 2myac which is quite natural for the composite 
cr-meson of the NJL model. However, we cannot prove equivalence of the two models when 
valence quarks are present. This is because the vacuum contribution to scalar density is 



strongly affected by the valence quarks and antiquarks. This effect is simply ignored in the 
sigma model. In contrast, the Dirac sea and valence quarks are treated self-consistently in 
the NJL model. In this respect it represents a more fundamental theory. 

Scaling expansion. -The evolution of a system out of equilibrium is governed by transport 
equations. For our purpose we use the relativistic Vlasov equation consistent with both the 
linear sigma model [|lOl and the NJL model 0, ^, |2^ : 



Pm 
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f{x,p) =0. 



(10) 



Here f{x,p) is the scalar part of the fermion distribution function. In general f{x,p) depends 



on 4-momentum p and 4-coordinate x. On the mass shell p^p^ 



m 



{x) and f{x^p) goes 



over into the sum of quark and antiquark occupation numbers as written in eqs. (^ and 
(1). The vanishing collision integral in the r.h.s. of eq. (p!0|) corresponds to a situation when 
collisions between quarks and antiquarks cease i.e. after thermal freeze-out. This might be 
a reasonable approximation for later stages of the expansion when the temperature drops 
below 140 MeV and the mean free path of quarks and antiquarks becomes long [|T6| . 

It is not trivial to find a solution of the Vlasov equation in a general case of (3 + 1) 
dimensional evolution. But in this letter we consider a simplified case of the spherical and 



homogeneous (Hubble-like) expansion [^. It is characterised by the collective 4- velocity 
y^fj. — x>^/t, where r = vi^—r^ is the proper time and r = {x,y, z). Assuming that m{x) 
depends on one variable r only, one can show [^ that a general solution of the Vlasov 
equation ([10|) is any function of the scaling variable t|p|. This function can be specified 
under the assumption that at the time of freeze-out, t = tq, the quark and anti-quark 
occupation numbers are given by the equilibrium Fermi-Dirac distributions. Then, at r > tq 
the quark and antiquark occupation numbers can be obtained as: 



nq[,T,p) 



^(^p)2 + m2(ro)-/io 



exp 



Tn 



+ 1 



nq[,T, Pj 



Un 



J, P, /Uo 



-f^o), (11) 



where m(ro), Tq = T(ro) and /io = Ai(To) are the constituent mass, temperature and chemical 
potential at r = tq. Below we set the chemical potential to zero, /iq = 0, so that Uq = n-q. 
The initial state for the scaling expansion is chosen at the respective critical temperatures 
for both models, i.e Tq = Tc. 

In the case of homogeneous spherical expansion the d' Alembertian is reduced to df^d^ — *• 
-^ + f ^- Now equations (^ and (|^) are ordinary differential equations which can be solved 



numerically by standard methods E^. In numerical calculations the initial values of fields 
and their first derivatives were selected within the range determined by the root-mean-square 
values of corresponding thermal fluctuations at r = Tq. At times t > tq the field fluctuations 
do not stay in thermal equilibrium but rather "propagate" according to equations of motion. 
In the quench scenario, introduced in ref. 0], one considers thermal equilibrium ini- 
tially and then lets the fields evolve out of equilibrium in the zero-temperature potential 
U{a,7f) (p. The fields roll down from the unstable local maximum of the Mexican Hat 
potential at a ~ 0, vr ^ towards the stable minimum at a = /^ and vf = and oscillate 
around these values. This quench scenario corresponds to a sudden removal of the heat- 
bath which initially provides thermal and dynamical equilibrium. As pointed out in |^, a 
more realistic scenario is when the potential changes gradually during the evolution and 
finally approaches U{a,7i). Gavin and Miiller [0 and recently Randrup fl^ considered a 
model where the fields evolve in an effective potential generated by mesonic quasi-particles 



treated in the Hartree approximation (hot anneahng). Our picture is somewhat similar, 
but instead of mesonic quasi-particles we consider quarks and antiquarks. However, there 
exists an important difference between [^ [1^ and the present approach. In refs. ||^, |14| the 



quasi-particles are assumed in thermal equilibrium throughout the whole evolution, while 
in our approach this assumption refers to the initial state only. At later times quarks and 
antiquarks follow the scaling solution of the Vlasov equation which differs from the thermal 
equilibrium distribution. 

Numerical results,- Fig 1 shows the effective potential for different proper times as pre- 
dicted by our model. One clearly sees that the potential evolves quite slowly and approaches 
the zero-temperature potential f/(cr, vf) only at about Stq. The asymmetry of the potential 
is due to the symmetry-breaking term Ha in the Lagrangian. Note that the effective po- 
tential depends on the initial conditions for the fields through the constituent quark mass 
appearing in ([TTD . The behaviour of the effective potential for the NJL model is qualitively 
similar. 

Now let us discuss numerical solutions of the equations of motion for scaling expansion. 
In the case of the sigma model we consider an initial temperature To = 132 MeV and initial 



proper time tq = 7 fm/c [Ty, ^ . At this temperature the root-mean-square fluctuations 
of fields and their derivatives are 38 MeV and 35 MeV/fm, respectively. Fig 2 displays 
the meson fields together with the constituent quark mass as functions of proper time. 
One observes strong oscillations of both the constituent quark mass and the fields. The 
complicated shape of the curves for the a-field and quark mass is due to the coupling of 
oscillations in different directions in isospin space. The constituent mass and the a-field 
first cross their equilibrium values already at about 1.3ro, although the effective potential 
at this time is still far away from the assymptotic shape. Notice the strong long-wavelength 
oscillations of the vr-field that may lead to rather big DCC domains ||^. 

It is instructive to compare three different scenarios: quench [Q, annealing and the 
collisionless expansion as discussed above. Fig. 3 shows the cr-field evolution for these 
three cases. The initial conditions are all around zero. In the quench scenario we see a 
large overshoot in the a-field of about 33 MeV compared to the equilibrium value of /^, 
then it oscillates strongly about the equilibrium value. This is a direct consequence of the 
zero-temperature potential used in this case. The fields just "glide" down from the top 
of the Mexican Hat to the minimum and then "climb" up again due to the large kinetic 
energy. In the annealing scenario we clearly observe the effect of the gradual change of the 
effective potential. The fields go "softer" towards equilibrium. In our scenario the effect of 
the changing effective potential (see Fig 1) is quite similar, except for stronger oscillations 
and slightly faster equilibration. 

In fig. 4 we compare the proper time evolution of the constituent quark mass calculated 
within the linear sigma-model and the NJL model. The initial state for the NJL model is 
chosen at Tq = 190 MeV corresponding to tq = 5 fm/c. Curve a) is calculated by combining 
eq. (|TT1) with the standard gap equation for the NJL model, Ga{x) = —1. Since in this 
lowest order approximation the gap equation does not contain any derivative terms, the 
constituent quark mass evolves smoothly to the vacuum value m^jac- One can say that it 
just follows the minimum of the evolving effective potential. In contrast, when eq. (0) is 
combined with the generalized gap equation (^ the character of solution changes drastically 
(curve (b)). First of all, one sees much faster rise of the constituent mass right after the 
beginning of the non-equilibrium stage. Second, the approach to the asymptotic value goes 
through strong oscillations, i.e. very similar to the sigma model (curve (c)). In the NJL case 
the oscillations have slightly longer period and smaller amplitude than in the sigma model. 
This is partly related to the difference in the critical temperatures and the corresponding 



initial times. As already mentioned above, this difference comes from the different treatment 
of the Dirac sea and valence quarks in these two models. 

Summary.- We applied two chiral models with quark degrees of freedom, i.e. the linear 
sigma model and the NJL model, to study dynamics of the chiral transition in expanding 
quark-antiquark plasma. A generalized gap equation with derivative terms was used for 
the NJL model. The calculations were carried out for a spherical homogeneous system and 
comparison of three different scenarios, quench, hot annealing and collisionless expansion, 
was made. Qualitively similar oscillatory behaviour of the constituent quark mass was found 
in all cases. We conclude that strong space-time fluctuations of the chiral fields are char- 
acteristic for a rapid chiral transition and they can be used for its experimental detection. 
In the future we are planning to study DCC domain formation and clustering of quarks as 
well as dissipative effects and soft pion production in the course of the chiral transition 
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figure captions 



Figure 1: Effective potential as a function of a- field (vr = 0) for different proper times 
(indicated in the figure). The zero-temperature potential f/(o", 0) is also shown. 



Figure 2: Meson fields and constituent quark mass in (units of /tt) as functions of proper 
time (in units of tq = 7 fm/c) calculated for the sigma model. The initial conditions at 
r = To are: $ = (38, 0, 45, 0) MeV and |f = (40, 35, 0, 0) MeV/fm. 



Figure 3: Proper time evolution of the cr-field calculated in the sigma model for three 
different scenarios. Dotted curve corresponds to the quench. The other curves are calcu- 
lated for scaling expansion with mesonic quasi-particles (solid line) and quark- ant iquark 
quasi-particles (bold line). The initial conditions are all zero 



Figure 4: Constituent quark mass as a function of proper time calculated within the NJL 
model without (a) and with (b) derivative terms. For comparison the result for the linear 
sigma-model is also shown (c). Explicit symmetry breaking terms are included and zero 
initial conditions are imposed in cases (b) and (c). 
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